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ABSTRACT 

Context. Observations and analysis of solar-type oscillations in red-giant stars is an emerging aspect of asteroseismic analysis with 
a number of open questions yet to be explored. Although stochastic oscillations have previously been detected in red giants from 
both radial velocity and photometric measurements, those data were either too short or had sampling that was not complete enough 
to perform a detailed data analysis of the variability. The quality and quantity of photometric data as provided by the CoRoT satellite 
is necessary to provide a breakthrough in observing p-mode oscillations in red giants. We have analyzed continuous photometric 
time-series of about 1 1 400 relatively faint stars obtained in the exofield of CoRoT during the first 150 days long-run campaign from 
May to October 2007. We find several hundred stars showing a clear power excess in a frequency and amplitude range expected for 
red-giant pulsators. In this paper we present first results on a sub-sample of these stars. 

Aims. Knowing reliable fundamental parameters like mass and radius is essential for detailed asteroseismic studies of red-giant stars. 
As the CoRoT exofield targets are relatively faint (11-16 mag) there are no (or only weak) constraints on the stars' location in the 
H-R diagram. We therefore aim to extract information about such fundamental parameters solely from the available time series. 
Methods. We model the convective background noise and the power excess hump due to pulsation with a global model fit and deduce 
reliable estimates for the stellar mass and radius from scaling relations for the frequency of maximum oscillation power and the 
characteristic frequency separation. 

Results. We provide a simple method to estimate stellar masses and radii for stars exhibiting solar-type oscillations. Our method is 
tested on a number of known solar-type pulsators. 

Key words, stars: oscillations - stars: fundamental parameters - techniques: photometric 



1. Introduction 

Stars cooler than the red border of the instability strip have con- 
vective envelopes where turbulent motions act over various time 
scales and velocities (up to the local speed of sound), producing 
acoustic noise which can stochastically drive (or damp) reso- 
nant, p-mode oscillations. All cool stars with convective outer 
layers potentially show these solar-type oscillations typically 
with small amplitudes. The oscillation amplitudes are believed 
to scale with the luminosity and are, therefore, more easily ob- 
served in evolved red giants than in main-sequence stars, which 
opens up a promising potential for asteroseismic investigations 
of evolved stars. Their larger radii, however, adjust the pulsa- 
tion periods from minutes to several hours to days. This in turn 
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complicates ground-based detection and calls for long and unin- 
terrupted observations from space. 

It is believed that the global characteristics of solar-type os- 
cillations, like the frequency range of pulsation or their ampli- 
tudes, are predetermined by the global properties of the star, like 
its mass or radius. It should therefore be possible to deduce the 
stellar fundamental parameters of a solar-type pulsator from the 
global properties of the observed os cillations. Recent investiga- 
tions in this context were made by Gilliland (2008) ,who an- 
alyzed pulsation amplitudes and timescales in several hundred 
red giants in t he galactic bulge o bserved by the Hubble Space 
Telescope and Stello et al. (2008), who determined asteroseis- 
mic masses for eleven bright red giants observed with the star 
tracker of the WIRE satellite. In this paper we measure global as- 
teroseismic quantities for 3 1 red giants observed with the CoRoT 
satellite and use well-known scaling relations to estimate their 
masses and radii. 

The satellite CoRoT (HaalirJ |2006l) is continuously col- 
lecting three-color photometry for thousands of relatively faint 
(about 11-16 mag) stars in the so-called exofield with the pri- 
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mary goal to detect planetary transits. But the data are also per- 
fectly suited for asteroseismic investigations. And indeed, a first 
processing of the C oRoT exofield data reveals a variety of os- 
cillating red giants. Hek ker et alj d2009l) report on the clear de- 
tection of solar-type oscillations in several hundreds of red gi- 
ants among the ~1 1 400 stars observed during the first 150 days 
CoRoT long-run (LRcOl) campaign. 

First results for oscillations i n red giants observed in the 
CoRoT exofield are presented in |Pe Ridder et alj (|2009), who 
also discuss the question of the existence of non-radial modes in 
red giants with moderate mode lifetime s versus the presenc e of 
short living radial modes only. Briefly, iBarban et alj d2007l) in- 
terpreted the signal found in the MOST observations of the red 
giant e Oph as radial modes with relatively broad profiles corre- 
sponding to short mode lifetimes of ~2.7 days. Their result was 
consistent w i th wh at was found for the similar red giant £ Hya 
dStello et all 120061) . On the other hand, iKallinger et al] 12008) 
re-examined the same data set and found radial and non-radial 
modes with significantly longer mode lifetimes (10-20 days), 
which s upports the lifetimes infer red from theoretic al considera- 
tions by Houdek & Gough ( 2002) and only recently Dupre fet alj 
(2009). Additionally, the frequencies are consistent with those of 
a red-giant model that matches e Ophs' position in the H-R di- 
agram. So the actua l interpretation of the av ailable observations 
was unclear. It was iDe Ridder et al.l ([2009) who first provided 
unambiguous evidence for regular p-mode patterns of radial and 
non-radial modes with long lifetimes. 

The low-degree, high-radial order p modes of solar-type pul- 
sators are approximately equ ally spaced and are believed to fol- 
low an asymptotic relation dTassoulj ll980i). For a detailed as- 
teroseismic analysis, such as the comparison of the individual 
observed frequencies and/or their separations with those of stel- 
lar models, is it often important to constrain the parameter space 
by using reliable fundamental parameters. Apart from pulsation, 
solar-type pulsators also show significant power from the turbu- 
lent fluctuations in their convective envelopes opening the pos- 
sibility to study convective time scales and amplitudes. Without 
knowing the position of the star in the H-R diagram such en- 
deavours cannot be properly realized. 

This knowledge is indeed lacking for most of the 
faint CoRoT exofield stars. The available broadband color- 
information allows at best, and only in some cases, a rough de- 
termination of the effective temperature, but is not suitable to 
distinguish giants from main-sequence stars. We therefore try to 
extract the fundamental parameters from the time series alone, 
which is the main topic of this paper. We model the convec- 
tive background noise and the power excess due to pulsation 
with a global fit, which allows us to measure the so-called fre- 
quency of maximum oscillation power. With this and the large 
frequency separation we derive the stellar mass and radius from 
well-known scaling relations. As a first estimate we also obtain 
effective temperatures and luminosities from a comparison with 
evolutionary tracks. 

2. Observations 

Although the CoRoT satellite and the exofie ld data are well 
explained by others dDe Ridder et all 120091 and references 
therein), we will briefly summarize the instrument and data sets. 

CoRoT houses four lk x lk pixels CCD photometers fed by 
a 27-cm afocal telescope. The satellite's low-Earth polar orbit 
(period ss 100 min; * 167 //Hz) enables uninterrupted observa- 
tions of stars in its continuous viewing zones (two cones with 
~10° radius centered on the galactic plane at right ascension of 




Fig. 1. Marginal distributions of the Gaussian parameters (Eq.|2]) 
for the red giant "A" as computed by the MCMC algorithm. 
Median values and lcr limits are indicated by vertical solid and 
dashed lines, respectively. 

about 6:50h and 18:50 h, respectively) for up to six months. A 
summary of the mission is given in the pre-launch proceedings 
of CoRoT published by ESA (SP-1306, 2006). The two core sci- 
ence objectives, asteroseismology across the H-R diagram and 
the detection of transiting extra-solar planets, are tracked simul- 
taneously with two of the four detectors each. In order to check 
the color independency of presumed planetary transits, the stel- 
lar light is dispersed by a prism before it reaches the exofield 
detectors. In this paper we concentrate on the white light flux 
measurements which are obtained by adding the flux of the three 
color channels to improve the photometric quality. 

During the first 150 day long-run campaign, CoRoT pointed 
toward the coordinates (a, 6) = (19. 4h , 0.46°) from May to 
October, 2007, and gathered time series of about 1 1 400 stars 
sampled with a cadence of either 512 s or 32 s, depending on 
the predefined status of the star (the limited downlink capac- 
ity does not allow to sample all stars with the short cadence). 
Typically, each time series consists of about 25 000 or rather 
400000 data points with a duty cycle of more than 90%. We 
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Fig. 2. Original (light-grey) and heavily smoothed (dark-grey) 
power density spectrum of the CoRoT exofield time series of the 
red giant A and a global model (black line) fitted to the power 
density spectrum. The model is a superposition of white noise 
(horizontal dashed line), three power law components (dashed 
lines) and a power excess hump approximated by a Gaussian 
function. The dotted line indicates the model fit plotted without 
the Gaussian component. 
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Fig. 3. Fourier amplitude spectra for a sample of red giants observed by CoRoT. Black lines indicate a global model fit, and dotted 
lines show the model plotted without the Gaussian component and serve as a model for the background signal. The center of the 
Gaussian is adopted to be the frequency of maximum oscillation power, v max - Error bars indicate the + lcr uncertainties and diamond 
symbols correspond to the weighted mean frequency after correcting for the background signal. Note the different amplitude scales 
in the different panels. Inserts compare the comb response functions (grey) and marginal distributions (black) used to determine the 
large frequency separations. 



use the N2 data format (Sama di et aTl l2007). which is the output 
of a standard data reduction procedure, and detect and remove 
occasional jumps in the time series (caused by high energy par- 
ticles) and apply an outlier correction. In a next step, we compute 
Fourier power spectra of all time series and extract parameters 
which we believed to be characteristic for red giant pulsators, 
like the 1/f 2 characteristic or the existence of a power excess 
hump. Based on a pre-selection with these parameters, we use 
a semi-automatic routine to identify the pulsating red giants. A 
more detail ed description of ho w to identify the red giant stars 
is given in iHekker et alj d2009h . The data are available at the 
CoRoT download pag43- 

3. Power spectra modeling 

The turbulent motions in the convective envelopes of cool stars 
act on a similar time scale as the acoustic oscillations and poten- 
tially complicate the detection and analysis of solar-type oscilla- 
tions. Although the convective signal is stochastic, it follows par- 
ticular characteristics. It can be shown that such quasi-stochastic 
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variations cause correlations of consecutive measurements with 
the strength of the correlations exponentially decreasing for in- 
creasing time-lags. The Fourier transform of this correlated col- 
ored "noise" follows a power law, characterized by an amplitude 
and a characteristic frequency (or inverse time scale). 

For the Sun, it is common practice to model the background 
signal with power laws to allow accurate measurements of the 
solar osci llation parame t ers. Power law mode ls were first intro- 
duced by Harvevl(fl 985). Aig rain et alj d2004l) and only recently 
iMichel et al.l d2009b use the sum of power laws: P(v) = 
with P, = a;£?T;/(l + (27!t,v) c to fit the solar background, with 
v being the frequency, t, the characteristic time scale, and C, the 
slope of the power law. a, serves as normalization factor for & = 
f Pi(v)dv, which corresponds to the variance of the stochastic 
variation in the time domain. The slope of the power laws was 
origi nally fixed to 2 in H arvey's models, but lAigrain et al.l(l2004l) 
and Michel et al. (2009) have shown that, at least for the Sun, the 
slope is closer to 4. The number of power law components usu- 
ally varies from two to five, depending on the frequency cover- 
age of the observations. Each power law component is believed 
to represent a different class of physical processes such as stel- 
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lar activity, activity of the ph otospheric/chromos ph eric magnetic 
netwo rk, or granulation (see Aig rain et alj|2004l or lMichel et alj 
2009 and references therein) with time scales for the Sun ranging 
from months for active regions to minutes for granulation. 

First tests with power law fits to the CoRoT photometry have 
shown that the presence of an additional power due to pulsations 
significantly distorts such a fit and requires an additional compo- 
nent to model the entire spectrum. Since the shape of the pulsa- 
tion power excess seems to be well approximated by a Gaussian, 
we model the observed power density spectra with a superposi- 
tion of white noise, the sum of power laws, and a power excess 
hump approximated by a Gaussian function 

P(v) = P„ + Y ^ + P. ■ e -C"--*) 2 /^), (l) 

where P„ represents the white noise contribution and A, and B, 
are the amplitudes of the stellar background components and 
their characteristic frequencies (or inverse time scales), respec- 
tively. The frequency coverage which results from the 150d 
CoRoT observations is sufficient to use three power law com- 
ponents. P g , v max , and cr g are the height, the central frequency, 
and the width of the power excess hump, respectively. 

We use a Bayesian Markov-Chain Monte Carlo (MCMC) al- 
gorithm to fit the global model to the power density spectrum. 
The algorithm samples a wide parameter space and delivers 
probability distribu tions for all relevant qua ntities. The proce- 
dure is described in Gruberba uer et al.l d2009l) and was originally 
designed to fit Lorentzian profiles to the p-mode spectrum of a 
solar-type pulsator. Our problem is quite similar, and it was triv- 
ial to adapt the code in order to fit a global model to the power 
density spectra instead of a sequence of mode profiles. The ad- 
vantage of the algorithm is its stability and insensitivity to wrong 
initial parameters, and also that it delivers reliable parameters as 
well as realistic uncertainties in a fully automatic way. For the 
frequency parameters (Z?,) we have sampled the entire frequency 
range of interest from to 150//Hz. For higher frequencies, the 
power spectra are potentially contaminated by instrumental arti- 
facts due to the satellites' orbital period. The amplitude parame- 
ters (P n , A,, and P g ) were allowed to vary from zero to the high- 
est amplitude peak in the spectrum. Only v max and cr g were kept 
within reasonable limits (0.5 to 2 times the value we inferred 
from a visual inspection of the spectra). After some 500000 it- 
erations we calculated the most probable value and its lcr uncer- 
tainty for all fitted parameters from their marginal distribution 
and constructed the most probable global model fit. 

As an example, we show in Fig.Q]the marginal distributions 
of the Gaussian parameters as computed by the MCMC algo- 
rithm for star A. The corresponding most probable global fit 
is given in Fig. [2] along with the original and heavily smoothed 
power density spectrum, the white noise, and power law com- 
ponents of the fit. Fig. [3] shows a sequence of amplitude spectra 
selected from the 3 1 analyzed red giants with the corresponding 
global model fits. Note that the fits are calculated in power but 
are presented in amplitude for better visibility. The sequence im- 
pressively demonstrates that the frequency and amplitude range 
of the oscillations scale with the time scales and amplitudes of 
the background signal. 

In a next step we use the white noise and power law com- 
ponents of the global model (doted lines in Fig.[2]and[3]l to cor- 
rect the power density spectra for the background signal, which 
should then include only the oscillation signal. But what is more 
interesting in this context, we assume the center of the Gaussian 
to locate the centroid of the power excess hump, and as long as 



the power excess hump is symmetric, the center of a Gaussian fit 
should equal the frequency of maximum oscillation power. To 
test this assumption we compute the weighted mean frequency, 
v, in the frequency range of pulsations (v max ±5cr g ), where we 
use the residual power after correcting for the background sig- 
nal as weight. Although for seven of the eight stars shown in 
Fig. [3] v is within v max ±lcr, v seems to be systematically shifted 
towards higher frequencies. Indeed, for 19 out of the 31 stars in 
our sample, v is higher than v max . But as the average shift for the 
19 stars (0.6%) is about four times smaller than the stochastic 
error of v max , we expect the systematic error on our estimate of 
Vmax to be negligible. Note that for the Su n the effect can not b e 
neglected. Based on SOHO/VIRGO data dFrohlich et al.Lll997l) . 
we find v to be shifted by almost 2% towards higher frequencies. 

The second interesting parameter that can directly be de- 
termined from the observed power spectrum is the large fre- 
quency separation, Av, of consecutive radial overtone modes of 
the same spherical degree. Since this frequency sep aration is, at 
least for the Sun (see, e.g., iBroomhall et all 120091) . a function 
of the frequency itself, it is not straightforward to identify an 
average value for all observed modes. Such a value depends on 
the actual number and frequency range of the observed modes 
and is difficult to compare for different observations. We there- 
fore specifically chose to define Av as the average frequency 
separation in the frequency range of the maximum oscillation 
power. To identify Av we use ag ain the Bayesian MCMC algo- 
rithm (Gruberbauer et all 120091) But instead of fitting a global 
model, we fit a sequence of equidistant Lorentzian profiles to the 
power density spectra 

_ 2 1 + 4[v - (v + Av ■ i/2)] 2 ■ (tit) 2 ' (2) 

where a,- is the rms mode amplitude of the z'-th profile, vo cor- 
responds to the frequency of the central mode, Av is the spac- 
ing, and t is the mode lifetime which is assumed to be equal 
for all five modes. In principle we do not expect all modes to 
have the same lifetime, but this is a marginal assumption which 
significantly stabilizes the fit and has very little impact on the de- 
termination of Av. The model obviously represents either three 
consecutive radial modes with two interjacent dipole modes or 
three consecutive dipole modes with interjacent radial modes, 
depending on which mode the sequence is centered. In either 
case, the fitted spacing corresponds to the large frequency sep- 
aration in the frequency range where the maximum oscillation 
power is seen. In our analysis the MCMC algorithm is allowed 
to vary the central mode frequency within v max +cr ? . Av and the 
mode lifetime are sampled between 0.5 /vHz and 2cr g and 1 and 
100 days, respectively. The individual mode amplitudes are al- 
lowed to vary between zero and four times the highest amplitude 
peak in the spectrum. After half a million iterations the algorithm 
delivers probability distributions for all fitted parameters and we 
calculate the most probable values and their uncertainties from 
their marginal distributions. 

The advantage of our method comp ared to, e.g., the comb- 
response function (Kiel dsen et all 1 19951) or an autocorrelation 
spectrum is that it takes the Lorentzian-like form of the signal 
into account and is therefore less sensitive to the stochastic na- 
ture of the signal. Examples for the residual power density spec- 
tra and the most probable fits are shown in Fig. [4] Interestingly, 
the presence of additional modes which are not taken into ac- 
count in our model does not influence the fit. This can be seen for 
instance from the power density spectrum of star A, where the 
MCMC algorithm correctly identifies the I - and 1 modes and 



Kallinger et al.: Oscillating red giants in the CoRoT exo-field 



5 



x 




w — 

C 15 - 



S 




Frequency (|iHz) 

Fig. 4. Residual power density spectra for star A and B after cor- 
recting for the background signal. The most probable model-fits 
used to determine Av are indicated by black lines. 



does not consider the additional peaks at about 72 and 79 yuHz, 
which are most likely I - 2 modes. We compare the marginal 
distributions for Av from our MCMC algorithm with the comb- 
response functions (both with arbitrary ordinates) in the inserts 
of Fig. [3] Although for some stars both methods give consistent 
results, the values can differ by more that 0.2 //Hz. 

We expect this ambiguity to be due to the stochastic nature 
of the signal. The observed time series represent a single real- 
isation of a damped and stochastically excited signal. Another 
realisation might result in a different measurement for Av and 
Vmax- In order to examine this we simulate different realisations 
of the same solar-type oscillation signal. We use the original 
time series of star A and pre-withen all significant peaks in the 
frequency range of pulsation. The residual time series should 
now include only the intrinsic background signal. We then gen- 
erate 250 time series following the procedure in IChaplin et al.l 
(119971) with each data set representing a different realisation of 
four radial orders of equidistant / = 0, 1, and 2 modes with Av = 
7. 20 //Hz and an arbitrary value for 6v - v n j - v„_i / + 2 = 1 //Hz. 
The time domain rms amplitudes of radial and / = 1 modes are 
set in a way that the mode heights follow the Gaussian shape 
of the original power excess hump. The rms amplitudes of / = 
2 modes are set arbitrarily to half the value of the closest radial 
mode. The largest amplitude mode is centered on 76.2//Hz, and 
all modes have a lifetime of 20 da ys, which is a t ypical value for 
intermediate luminous red giants (Dupret et all 120091) and cor- 
responds to what we have found for the intrinsic modes. Finally 
we superpose each simulated data set with the residual time se- 
ries, calculate the power density spectrum, and apply our algo- 
rithms to determine v max and Av. The simulations give an av- 
erage value v max = 74.42 //Hz with a rms scatter of 0.55 //Hz, 
which is well within the 1 <x uncertainty of the originally deter- 
mined value (74.32+0.81 //Hz). The situation is different for Av. 
If we determine Av from the frequency of the largest peak of the 
comb-response function, then the stochastic nature of the signal 
adds a rms scatter of about 0.21 //Hz in our simulation. This is 
much compared to the rms scatter of only 0.04 //Hz if we use our 
Bayesian MCMC approach, which is again compatible with the 
lcr uncertainty of the original value (0.04 //Hz). This assures us 



that the stochastic nature of the oscillation signal adds no signif- 
icant additional uncertainty in our subsequent analysis. We list 
v max , Av and the corresponding lcr uncertainties for our sample 
of red giants in Table|2] 

4. Asteroseismic determination of fundamental 
parameters 

Stellar masses of field stars are usually determined by compar- 
ing the location in the H-R diagram with evolutionary tracks. But 
when stars evolve to become red giants, their evolutionary tracks 
move together and a relatively narrow range in the H-R dia- 
gram covers a large range in mass. Consequently, the mass deter- 
mined from a location in the H-R diagram becomes quite uncer- 
tain. Additionally, the actual position and slope of the red-giant 
branch in stellar evolutionary calculations depend very much on 
the parameters used to compute the models. It is, however, be- 
lieved that the global properties of solar-type oscillations, like 
the frequency range where they can be observed, depend on the 
fundamental parameters of the star. It should therefore be possi- 
ble to determine these parameters from global properties of the 
observed oscillations. 

The amplitudes of solar oscillations are modulated by a 
broad envelope with its maximum at a frequency of about 3 mHz. 
The center and the shape of the envelope is defined by the ex- 
citation and damping where the later can b e assu med to be 
G aussian (see previous sectio n). iBrown et"ail (Il99ll) and latter 
on Kieldsen & Bedding ( 1995b have shown that the frequency of 
maximum oscillation power, v max , of p-mode oscillations scales 
to good approximation with the acoustic cutoff frequency, which 
sets limits on the maximum frequency for acoustic oscillations. 
They predict v max by scaling from the Sun as 

v raax = (M/MQ) ■ (R/Roy 2 ■ (r eff /5777 Ky l/2 ■ 3050//Hz. (3) 

It has been shown that this simple scaling relation gives very 
good estimates for the freque ncy of maximum oscillation power 
for less evolved stars (e.g. B edding & Kieldsenl 120031) . but it 
cannot a pr iori be assumed th at it holds also for stars of the gi- 
ant branch. Stello et al.l d2008l) . however, have demonstrated for 
a number of bright red giants observed by the star tracker of the 
WIRE satellite that the measured v max is in reasonable agreement 
with the values that can be expected from their fundamental pa- 
rameters. In their analysis they use the luminosities and effective 
temperatures determined from Hipparcos parallaxes and infrared 
photometry, respectively, to substitute the radius in Eq.[3]accord- 
ing to L oc R 2 ■ T 4 . 

Our sample of red giants is far too faint to measure paral- 
laxes. But the CoRoT observations are significantly better than 
the WIRE observations in terms of duration and precision, which 
enables us to also extract the larg e frequency separation. For 
the Sun, Toutain & Frohlich (1992) have determined a large fre- 
quency separation of about 134.92 //Hz at the radial order where 
the maximum oscillation power is seen (n = 21). The large 
frequency separation reflects essentially the global properties 
of the star and is believed to scale with the dynamical time 
scale and therefore with th e square root of the mean density. 
Kieldsen & Bedding] d 19951) predicted Av by scaling from the 
Sun as 

Av = (M/MQ) 112 ■ (R/RQ)- y2 ■ 134.92 //Hz. (4) 

Knowing v max and Av (and T e fr), it is now easy to derive the 
stellar mass and radius: 

R/RQ = (v,nax/v max , ) ■ (Av/Avor 2 • (7W/5777 K) 1 ' 2 (5) 
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Table 1. Mass and radius of stars used to test our asteroseismic mass and radius determination approach. Mass and radius as taken 
from the literature are given in brackets. 
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Av 
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(R 
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IjuHz] 






solar units 




Arcturus 


3.47±0.03 


0.825+0.05 


4290+30 


27.9+3.4 


0.8+0.2 


(25.4+0.3 


0.7-2.0) 


HD181907 


29.1+0.6 


3.47+0.12 


4760+65 


13.1+1.0 


1.5+0.23 


(12.3+0.6 


1.5-2.0) 


jSOph 


46.0±2.5 


4.1+0.2 


4470+100 


14.9+1.7 


3.1+0.8 


(12.2+0.8 


<3.2) 


eOph 


53.5±2 


5.2+0.1 


4877+100 


10.7+0.6 


1.8+0.3 


(10.4+0.5 


2.02) 


f Hya 


92.3±3 


7.0+0.2 


5010+50 


10.2+0.7 


2.9+0.4 


(10.4+0.5 


2.93-3.15) 


M67 13 


208.9±4 


15.9+0.2 


4966+50 


4.6+0.1 


1.33+0.1 


(4.3+0.3 


1.35) 


v Ind 


313+10 


24.25+0.25 


5300+100 


3.04+0.12 


0.91+0.1 


(2.96+0.13 


0.81-0.89) 


f Her A 


700±50 


43+0.5 


5825+50 


2.3+0.17 


1.2+0.26 


(2.51+0.11 


1.3-1.5) 


/3Hyi 


1020±50 


56.2+0.2 


5872+44 


1.94+0.1 


1.28+0.2 


(1.81+0.02 


1.0-1.2) 


HD49933 


1657±28 


85.2+0.5 


6500+75 


1.44+0.03 


1.21+0.07 


(1.46+0.05 


1.325) 


/j Ara 


1900+50 


90+1 


5813+40 


1.40+0.05 


1.23+0.1 


(1.36+0.03 


1.1-1.3) 


a Cen A 


2410+130 


106+0.2 


5770+50 


1.28+0.07 


1.30+0.2 


(1.26+0.03 


1.124+0.008) 


aCenB 


4090+170 


161.4+0.06 


5300+50 


0.90+0.04 


1.04+0.13 


(0.85+0.02 


0.934+0.007) 



M/MO = (R/RO) 3 • (Av/Av©) 2 



(6) 



We have tested this method for a number of well-known solar- 
type pulsators and compare in Table[TJthe seismic masses and 
radii with independent measurements given in the literature. The 
specific values for the latter are taken from 

- Arcturus: Vm ax and Av are given bv lTarrant et alj (12007!) and 
iRetter etaTI d2003l) , respectively. The ra dius is based on the 
Hipparcos para llax ( van Leeuwenll2007l) and interferometric 
measurements (lLacour et al. , 2008). The mass range is esti- 
mated from the average surface gravity listed in the VizieRQ 
database (logg = 1.72+0.2) and the interferometric radius. 

- HD181907: We determine y max from the CoRoT observa- 
tions published in Carri er et al.l d2009l) . who also derived Av. 
The mass is estimated from a comparison between the star's 
position in the H-R diagram and metal-poor evolutionary 
tracks. 

- p Oph: v max and Av were determined from unpublished 
MOST ph otometry. The radi us is based on the Hipparcos 
parallax (van Leeuwen, 2007) and interferometric measure- 
ments ( Richic hi et al.L 12005). An upper limit for the mass 
is estimated from the average surface gravity (logg = 
2.42+0.3) taken from the VizieR database and the interfer- 
ometric radius. 

- e Oph: v max and Av are given bv lKallinger etaT] d2008l) . who 
determined the mass from a detailed comparison of observed 
and mod el frequencies. The ra dius is based on the Hipparcos 
paralla x dvan Le euwen. 2007) and interferometric measure- 
ments dRichichi et al.Ll2005l) . 

- £ Hya : v m ax> Av, and the mass are taken from Fra ndsen et alj 
(2002) , where they estimate the mass from a compari- 
son between the star's position in the H-R diagram and 
solar-calibrated evolutionary tracks. We derive v max from a 
weighted average of the published frequencies. 

- M67 13: v max and Av are determined from the photometric 
time series kindly provided by D. Stello. The mass is esti- 
mated fronn^ochronefits to the color-magnitude diagram of 
M67 (IStello et alll2007h. 



v Ind: Av is given bvlCarrier et a i] d2007l) . y m ,v is taken from 
iBedding et alj d2006l) . who also provide a mass range based 
on a comparison between the star's position in the H-R dia- 
gram and evolutionary tracks. 



- Z" Her A: v max , Av, and the mass are taken from lMartic et al.l 

Eoolh . 

- fi Hyi: v max and Av a re given bv lKieldseneTail d2005l) and 
Bedd ing et al.l (1200 ll) . respectively. The radius is bas ed on 
interferometric measurements from lNorth et alj (|2007), who 
also provide a summary of non-seismically determined val- 
ues for the mass. 

- HD49933: v max and Av are given bv lKaTlinger et al.l(l2009bl) . 
who determined the mass from a detailed comparison of ob- 
served and model frequencies. 

- n Ara : v max and Av were extracted from iBouchv et al.l 
d2005l). A non-seism ic estimate for the mass can be found 



Bazot etall (120051) . 
a Cen A & B: v max and Av are tak en from Kieldse n et al.l 



( 2005 ). Masses were determined by Guent her & Dem arque 
(2000) using Hipparcos parallaxes and the binary mass ratio. 



http://vizier.u-strasbg.fr/viz-bin/VizieR 



If not explicitly mentioned the radius is determined accord- 
ing to L oc R 2 ■ rf ff with the effective temperature and luminosity 
taken from the listed reference (or references therein). Although 
v max is in some cases only a rough estimate from a published 
power spectrum, our approach yields quite accurate masses and 
radii in a large portion of the H-R diagram. For faint stars it often 
turns out that no or only poor estimates for the effective temper- 
ature are available. This is, however, not very critical for stars 
on the giant branch, as v max depends only on the square root of 
Teff, and pulsating red giants are expected to populate only a rel- 
atively narrow temperature range (~4200 to 5300 K). It should 
therefore be possible to get a reasonable asteroseismic mass and 
radius for red giants even if an accurate T e ff is not available. This 
is shown for the red giants in TablefTJ (Arcturus to £ Hya), where 
we ignore the known values for T^g and fix the temperature to a 
typical value of 4750 K. 

For our CoRoT sample of red giants we first esti- 
mate the effective temperature from 2MA SS photomet r ic col- 
ors and color-temperature calibrations (|Masana et al., 2006; 
Gonzalez Hernandez & Bonifacio, 120091) . Unfortunately, both 
calibrations result in temperatures which are syst ematically too 
cool. T he average temperature resulting from the iMasana et al.l 
(2006) calibration is about 3820 K, which would make our sam- 
ple of red giants either extremely metal rich or would put all 
the stars high up on the giant branch. Both explanations are 
not very plausible. We expect the discrepancy to be due to 
severe reddening which is difficult to estimate for such faint 
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Fig. 5. Theoretical H-R diagram showing the location of the stars used to test our asteroseismic mass and radius determination 
approach. Grey-filled dots (total sample) and black dots (stars presented in Fig.0 indicate the analyzed CoRoT pulsating red giants 
where the actual position in the HR-diagram is based on a comparison with evolutionary tracks. The errors bars correspond to 
the approximate uncertainties of our method and are significantly larger than the observationally based errors. Dashed black lines 
indicate isopleths for Av = 2.5, 5, and lO^uHz and dotted black lines show isopleths for v max = 20, 40, 80, and 160 ^Hz where we 
linearly interpolate for points between the evolutionary tracks. 



(an d therefore distant) stars. The situation i s sligthly better for 
the lGonzalez Hernandez & Bonifaciol (f2009) calibration, but we 
decided to ignore the effective temperatures determined from 
color- temperature calibrations. 

Instead, we use a different approach. We calculate an initial 
guess for the mass and radius from Eqs.[5]and[6]by fixing the 
effective temperature to a typical value of 4750 K. In a next step 
we compare the initial mass and radius to those of a grid of solar- 
calibrated red giant models. Interpolation in the grid gives a bet- 
ter estimate for the temperature, which is used as a new input for 
Eq.[5] After about three iterations the procedure converges to a 
certain location in the H-R diagram where the final locus in the 
H-R diagram is independent from the starting value for T e g as 
long as the initial value is kept within the temperature range of 
our model grid (-4000 - 5500 K). 

The resulting fundamental parameters and their uncertainties 
are given in Table|2] The errors are based only on the uncertain- 
ties for v max and Av and range from 1.5 to 8.2% and 3.9 to 17% 
for the radius and mass, respectively. 

The red-giant models used to estimate the effective tem- 
peratures fall along evolutionary trac ks computed w i th the 
Yale Stellar Evolutio n Code YREC dGuenfher et all 119921; 
iDemarque et al.ll2007l) . The evolutionary tracks were computed 
for an initial helium and metal mass fraction (Y, Z) = (0.28, 0.02) 
with the mixing-length parameter a = 1.8 set to approximately 
meet the Sun's position in the H-R diagram with a one-solar- 



mass model at roughly the solar age. Note, although the models 
are not exactly calibrated to the Sun, we refer to them as solar- 
calibrated in the following. A more detailed descript i on of the 
used model phys ics can be found in iKallinger et alj (120081) or 
iKallingerl (l2009al) and refer ences therein. 

In Fig. [5] we show our sample of red giants in the H-R dia- 
gram along with the test stars (Table[TJ and contours of constant 
Vmax and Av. We whish to emphasize that the actual locations of 
the red giants are specific to the model grid used to estimate their 
effective temperatures. We assume that all red giants are compa- 
rable to the Sun in terms of their initial chemical composition 
and mixing-length parameter. This might be correct for some of 
them but not for others. The effective temperature and luminos- 
ity of a model with a given mass and radius do not only depend 
on the model's age but also for instance on the initial chemical 
composition of the model. Metal-poor models, e.g., are shifted 
towards higher surface temperatures and luminosities compared 
to solar-abundant models. On the other hand, low and interme- 
diate mass red giants contract rapidly after they have reached 
the tip of the red giant branch (RGB) to settle on the He-core 
burning main sequence at a somewhat higher temperature be- 
fore they start to climb the asymptotic giant branch (ABG). In 
other words, a red giant with a given mass and radius is located 
at different positions in the H-R diagram depending on, e.g., the 
chemical composition and/or the evolutionary stage. 
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Table 2. Summary for the analysed sample of red gaints. V magnitudes are taken from the EXODAT database dDeleuil et aU l2006). 
Vmax an d Av are obtained from the photometric time series. Fundamental parameters are based on scaling relations for v max and Av 
and interpolation in a grid of solar-calibrated RGB models. The errors are based only on observational uncertainties. 





CoRoT-ID 


2MASS 


V [mag] 


Vmax [/"Hz] 


Av [pHz] 


R/R0 


M/MO 


(T eff 


L/LO) 


A 


101113062 


19263135-0002147 


13.46 


74.32±0.81 


7.20= 


=0.04 


7.69±0.11 


1.30±0.05 


(4670±10 


25: 


fcl) 


B 


101251252 


19272346+0119408 


13.72 


67.74±0.92 


6.50= 


=0.03 


8.62±0.15 


1.49±0.07 


(4684±14 


32: 


fcl) 


C 


101362522 


19281037+0049531 


13.90 


53.51±1.28 


5.53= 


=0.04 


9.32±0.26 


1.36±0.11 


(4612±21 


35: 


t3) 


D 


101034881 


19255782+0056022 


13.50 


45.59+0.90 


4.79= 


=0.02 


10.58±0.23 


1.49±0.09 


(4598+18 


45: 


t3) 


E 


101197556 


19270316-0032250 


13.35 


34.60±0.69 


3.79= 


=0.04 


12.75±0.40 


1.64±0.13 


(4558±23 


63: 


±5) 


F 


100838545 


19245182+0129416 


12.76 


31.03±0.65 


4.15= 


=0.05 


9.38±0.29 


0.78±0.06 


(4395±16 


29: 


±2) 


G 


101649216 


19301425-0007224 


12.46 


19.30±0.37 


2.82= 


=0.11 


12.48±1.00 


0.85±0.14 


(43 11 ±26 


48: 


t9) 


H 


101378942 


19281895+0105256 


12.76 


14.03±0.44 


2.05= 


=0.03 


17.25±0.69 


1.18±0.13 


(4294±22 


91: 


±9) 




101290847 


19273800+0112570 


13.20 


61.13±0.60 


5.93= 


=0.11 


9.33+0.35 


1.57+0.12 


(4676±21 


37: 


t3) 




101509360 


19291255-0011005 


13.06 


56.77±0.66 


5.74±0.03 


9.22±0.16 


1.42±0.06 


(4634±12 


35: 


±2) 




101081290 


19261943-0032595 


12.73 


51.59±0.42 


5.32= 


=0.06 


9.72±0.20 


1.43±0.07 


(4615+11 


38: 


±2) 




101041814 


19260095+0131378 


13.51 


43.40±0.57 


4.90= 


=0.06 


9.55±0.24 


1.15+0.07 


(4539±12 


35: 


t2) 




100886873 


19250710+0106045 


13.57 


40.89+0.75 


4.73= 


=0.04 


9.65±0.24 


1.10+0.07 


(4515+15 


35: 


t2) 




100483847 


19224577+0131126 


12.70 


39.15±1.07 


4.44±0.03 


10.50+0.30 


1.25±0.10 


(4534±19 


42±3) 




101232297 


19271621-0015395 


12.48 


40.02±0.91 


4.78= 


=0.03 


9.22±0.24 


0.98±0.07 


(4493±19 


31: 


±2) 




100974118 


19253501+0022085 


13.74 


39.03±0.66 


4.91= 


=0.06 


8.48±0.25 


0.81±0.06 


(4446±15 


25: 


±2) 




100716817 


19241108+0138078 


13.86 


37.40±0.77 


4.02= 


=0.03 


12.28±0.31 


1.65±0.11 


(4577+22 


59: 


±4) 




101218811 


19271113+0152095 


13.46 


34.53+1.10 


3.99= 


=0.04 


11.42±0.44 


1.30±0.14 


(4512+25 


49: 


±5) 




101242228 


19271998+0045148 


13.10 


31.17±0.70 


3.92= 


=0.08 


10.62±0.47 


1.01+0.10 


(4445±23 


40: 


±4) 




101136306 


19264008+0127108 


12.61 


32.22±0.58 


3.91= 


=0.04 


11.06±0.31 


1.14+0.08 


(4475±15 


44±3) 




101058180 


19260841+0151442 


12.61 


31.84±0.77 


4.08= 


=0.05 


9.98+0.35 


0.91±0.08 


(4430+21 


34: 


±3) 




100908597 


19251396+0034026 


13.95 


31.76±1.10 


3.93= 


=0.04 


10.78±0.43 


1.06±0.12 


(4459+26 


41: 


±4) 




101262795 


19272783+0100075 


14.16 


30.75±0.81 


3.92= 


=0.07 


10.49±0.47 


0.97±0.10 


(4435±24 


38: 


±4) 




101044584 


19260224-0019015 


13.00 


28.53±0.98 


3.95= 


=0.08 


9.49±0.51 


0.73±0.09 


(4368±24 


29: 


±4) 




101513442 


19291426+0009508 


14.32 


30.06+0.80 


3.64= 


=0.05 


11.91+0.46 


1.23±0.12 


(4473±21 


51+5) 




100855073 


19245704+0130376 


14.13 


26.37±0.52 


3.28= 


=0.02 


12.87±0.32 


1.26±0.08 


(4446±14 


58: 


±4) 




101150795 


19264557+0051573 


14.11 


25.05±0.61 


3.19= 


=0.07 


12.88+0.69 


1.19+0.15 


(4428±23 


57: 


±7) 




101654204 


19301756-0014319 


13.59 


26.71±0.90 


3.63= 


=0.03 


10.54±0.40 


0.85±0.09 


(4378±25 


37: 


±4) 




100998571 


19254285-0012038 


13.70 


24.29±0.52 


2.89= 


=0.04 


15.27±0.56 


1.64±0.14 


(4475±23 


84: 


±8) 




101449976 


19284824-0007431 


12.38 


21.50±0.56 


2.79= 


=0.02 


14.38±0.43 


1.27+0.11 


(4401±17 


70: 


±5) 




101029979 


19255560+0122554 


14.64 


21.85±1.05 


2.60±0.03 


16.95±0.89 


1.81±0.27 


(4474±44 


100± 


15) 



To illustrate this ambiguity we compare stellar models for a 
given mass and radius, but with different initial chemical com- 
position and during different evolutionary stages. The result is 
shown in Fig. [6] where we illustrate that 2.5 Mo RGB models 
with 10 Ro with a initial chemical composition of (Y, Z) = (0.25, 
0.01) and (0.32, 0.04) are about 155K hotter and 13% more lu- 
minous and about 135K cooler and 11% less luminous, respec- 
tively, than a solar-calibrated RGB model with the same mass 
and radius. Similar results can be expected for different mixing- 
length parameters. Whereas the parameterization of convection 
has only small effects on the surface properties of a star dur- 
ing early evolution, different mixing-length parameters result in 
quite different evolutionary tracks when the star ascends the gi- 
ant branch. This is because the mixing-length parameter sets the 
temperature gradient in the convective regions and thus con- 
trols how efficiently energy can be deduced from the interior. 
Consequently, the mixing-length parameter defines at what stage 
the star starts to climb the giant branch during the hydrogen shell 
burning phase. A slightly smaller effect can be expected for the 
different evolutionary stages. Our 2.5 Mo and 10 R© is about 
60K hotter and 5% more luminous in its AGB phase than the 
corresponding RGB model. Although most of our sample red 
giants are expected to have a mass below 2.5 Mo, we use here 
the more massive models. This is because YREC, just like other 
stellar evolution codes, is not able to follow the explosive He 
flash in low mass stars, and we currently have no AGB start- 



ing models available. We can follow the He core ignition only 
for higher mass models and evolve models on the AGB. We do 
not expect the effect to be significantly different for lower mass 
models, however. 

We have demonstrated in Fig. [6] that an unknown chemical 
composition, mixing-length parameter, and evolutionary stage 
adds a certain amount of uncertainty when we determine a star's 
position in the H-R diagram from its mass and radius. This am- 
biguity will add additional uncertainties on the asteroseismic 
masses and radii. To quantify this effect we have computed addi- 
tional sets of evolutionary tracks. For the sake of simplicity and 
as we expect the uncertainty to be smaller for an unknown evolu- 
tionary stage than for an unknown initial chemical composition 
we concentrate here on models with different initial chemical 
compositions, namely (Y, Z) = (0.25, 0.01) and (0.32, 0.04). In 
Table|3]we compare the fundamental parameters as they result 
from using the different model grids to estimate the effective 
temperature for Eq.|5] We list here only three stars which are, 
however, selected to cover the range of interest in the H-R di- 
agram. The spread in chemical composition adds an additional 
uncertainty of not more than about 1.8 and 5.5% on the radius 
and mass, respectively, which is smaller or at most comparable to 
the observationally based errors. The situation is different for the 
effective temperature and luminosity where the additional un- 
certainties are significantly larger than the observational errors. 
To account for this we have extended the error bars in Fig. [5] to 
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Fig. 6. YREC evolutionary tracks for 2.5 Mo models with dif- 
ferent chemical compositions (Y, Z) = (0.25, 0.01), (0.28, 0.02), 
and (0.32, 0.04) and mixing-length parameters a = 1.6, 1.8, and 
2.0. For the solar-calibrated track also the AGB part of the track 
is shown. The large dots indicate models with 10RO. 



Table 3. Fundamental parameters for the stars A, E, and H as 
they follow from using a metal-poor (MP), solar-calibrated (SC), 
or metal-rich (MR) grid to estimate the effective temperature. 
Also given are standard deviations in absolute units and percent. 



ID 


Grid 


R/RO 


M/MO 


Teffffl 


L/LO 


A 


MP 
SC 
MR 


7.82±0.12 
7.69±0.11 
7.57±0.12 


1.36±0.05 
1.30±0.05 
1.24±0.05 


4835+11 
4670+10 
4528+10 


30+1 
25+1 
22+1 




(T 


0.12(1.6%) 


0.06(4.6%) 


153(3.3%) 


4(15%) 


E 


MP 
SC 
MR 


12.98±0.39 
12.75±0.40 
12.56±0.38 


1.73+0.13 
1.64±0.13 
1.56+0.12 


4725+25 
4558+23 
4417+19 


75+6 
63+5 
54+4 




a 


0.21(1.6%) 


0.09(5.5%) 


154(3.3%) 


11(17%) 


H 


MP 
SC 
MR 


17.59±0.71 
17.25±0.69 
16.98±0.71 


1.25+0.13 
1.18+0.13 
1.13+0.12 


4462+46 
4294+22 
4160+28 


110+10 
91+9 
78+9 




a 


0.31(1.8%) 


0.06(5.0%) 


151(3.5%) 


16(17%) 



+ 150K and +17% for the effective temperature and luminosity, 
respectively. Note that these are still quite accurate fundamental 
parameters for such faint stars. 

From Eqs.[3]and|4]and the isopleths for v max and Av in Fig. [5] 
it is obvious that the two parameters are correlated to some ex- 
tent. Both parame ters depend on the stellar mass and radius. 
Stell o et all (12009) indeed found a tight relation between v max 
and Av fo r main-sequence an d red-giant stars, which was con- 
firmed by Hek ker et alj (120091) for CoRoT red giants. Although 
the relation seems to be very tight, no exact relation can be ex- 
plained from a theoretical point of view. It even turns out that 
the relation between v max and Av strongly depends on the stellar 
mass. This can be seen from Fig.[7j where we compare the ratio 
between Av and v max as a function of v max for models of different 
evolutionary tracks with the sample of red giants. The ordinate 
basically represents the i nverse radial order where the maximum 
oscillation power is seen ( Kieldsen & BeddingLll995h . We intro- 
duce this diagram as a sort of diagnostic asteroseismic diagram 
similar to a diagram which is usually used to estimate the mass 
and central hydrogen abundance for solar-type pulsators close 
to or on the main sequence fro m the measured large and smal l 
frequency separations (see e.g., iRoxburgh & Vorontsovl 120031) . 
Our diagnostic asteroseismic diagram is particularly useful on 
the giant branch and relies on an observational quantity, namely 
Vmax, which is easier to determine than the small frequency sep- 
aration. Other parameters like the helium core abundance can 
easily be added to this diagram. The diagram also demonstrates 
the relative robustness of our method to determine an asteroseis- 
mic mass and radius. Models of a given mass and radius (large 
dots) are only slightly shifted in the diagram if e.g. their initial 
chemical composition is significantly changed. 

We have illustrated to some extent what uncertainty can be 
expected for a red giant's position in the H-R diagram if the 
star's initial chemical composition and/or mixing-length param- 
eter and/or evolutionary stage is unknown. But there are also 
other effects in stellar evolution which carry a similar type of 
uncertainty. Examples are the overshoot parameter or a better 
description of convection than the MET which change the L-M- 



R-T e ff relations. These effects are difficult to estimate and are not 
on the scope of the current analysis. 

5. Conclusions and prospects 

We have shown that global properties of solar-type pulsations 
can be used to derive estimates for the stellar mass and radius 
by employing well-established and often used scaling relations. 
We have tested this approach on various prominent solar-type 
pulsators and applied it to a first sample of red giant pulsators 
observed by CoRoT. Despite the mentioned approximations the 
derived fundamental parameters can serve to constrain the start- 
ing values for a more detailed analysis. 

We note that we do not stop at this point. In a next step we 
will use the integral of the Gaussian part of our global power 
density model to deduce the total spectral power of solar-type 
pulsations, which we expect to scale with the luminosity-mass 
ratio. This, however, needs extensive calibration for solar-type 
pulsators with independently determined fundamental parame- 
ters, which we are currently carrying out. We believe that we 
can, for the first time, derive all basic fundamental parameters 
(mass, radius, luminosity, and consequently also the effective 
temperature) of a solar-type pulsatorby simply measuring global 
properties of its oscillations in the power density spectrum. This 
will have influence on various astrophysical applications. One 
can use it as a distance indicator, or one can study the behavior 
of convective time scales of stars as a function of their position 
in the H-R diagram. By comparing the individual pulsation fre- 
quencies with theoretical eigenfrequencies it should be possible 
to investigate the parameterization of convective models (e.g, the 
mixing-length parameter) in a region of the H-R diagram where 
stars are very sensitive to these parameters. 

Finally, we want to mention that we will apply our aster- 
oseismic fundamental parameter determination to all pulsating 
red giants observed by CoRoT, and also we plan to arrange an 
online database for them. 
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